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Abstract: We study parameters of an extended standard model. The neutrino sector 
is enlarged by one or two right-handed singlet fields and the Higgs sector contains one 
additional doublet. One-loop radiative corrections generate the mass for the light neutrino 
fields. The numerical analysis is performed varying the masses of heavy neutrinos and of 
the additional neutral Higgses. The parameters of the neutrino sector, allowing for the 
seesaw type-I mechanism, are restricted by experimental neutrino oscillation data. Both 
normal and inverted hierarchies of the light neutrino masses are discussed. 
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1 Introduction 

The precise interpretation of the neutral lepton fields in the particle physics Lagrangian 
is not settled yet, owing to the very small mass of the known neutrinos and the weakness 
of their interaction with other particles [1]. The observed neutrino oscillations support 
the notion that neutrinos have non-vanishing masses, calling for a modification of the 
Standard Model (SM). The size of the neutrino mass is not the only puzzle to solve. 
Absence of an electrical charge allows neutrinos to be their own antiparticles. The nature 
of the neutrinos - whether they are Dirac or Majorana particles - might be determined by 
future experiments. 

The Standard Model considers neutrinos as massless. Adding heavy right-handed 
neutral singlets and additional Higgs doublets, the authors of ref. [2] combined the seesaw 
mechanism (type-I) with the radiative mass generation. The spontaneous symmetry break¬ 
ing (SSB) of the SM gauge group leads to a Dirac mass term for neutrinos. The assumption 
that neutrinos are Majorana particles allows an additional term in the Lagrangian, namely, 
the Majorana mass term for the heavy singlets. 
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The model parameters allow small masses of the light neutrinos that are compatible 
with the experimental observations. We use this model in the formulation of Grimus and 
Lavoura [3, 4], restricting the number of additional Higgs doublets to one. We consider 
only 1-loop corrections to the neutrino mass matrix. The case of three additional heavy 
neutrinos added to three light neutrinos was studied e.g. in ref. [5, 6]. We assume either one 
or two heavy neutrinos. Our preliminary results were presented at several conferences [7- 
11]. Here we provide a more complete description of the performed numerical analysis. 

Our extended model has several subsets of the parameters. The neutrino sector is 
characterized by the masses of the heavy neutrinos (either one or two), and the strength 
of the coupling to the neutral Higgs fields. The masses of three light neutrinos are the 
result of our model parameters. They are subject to experimental constraints, namely 
the experimental neutrino mass differences, Attiq and as well as the experimental 

neutrino oscillation angles 612 , ^13, and 023 [12]. To estimate the neutrino oscillation angles 
from the neutrino mixing matrix we follow the ideas of ref. [13] . More details are given in 
appendix B. It should be noted that experimental data is usually interpreted in the “3 x 3” 
neutrino mixing model [1, 12], i.e. three flavoured neutrinos are considered as mixed states 
of three neutrino mass-eigenstates. We did not attempt to reinterpret the results in the 
context of an extended neutrino model. We expect the effect to be negligible. 

We parametrize the Higgs sector along the analysis of [14]. The Yukawa couplings 
are parametrized similarly to Grimus and Lavoura [3, 4], which coincide with [14] in the 
Higgs sector. For the numerical analysis we take the mass of the SM Higgs boson as 
rrifjo = 125 GeV [15] and allow the masses of two heavier Higgs bosons to vary in the range 
from 126 to 3000 GeV. 

The outline of the paper is following. Section 2 reviews the seesaw mechanism and 
the formalism of the two-Higgs-doublet model as used in our analysis. Sections 3 and 
4 describe our main results, namely, the calculated light neutrino mass spectra and the 
analysis of free parameters. Our findings are summarized in section 5. For completeness, 
the appendix section A describes the features of the weight vectors hi that relate the scalar 
Higgs fields to their mass eigenfields, and section B gives the details of the oscillation angle 
calculation. 


2 Description of the model 

We discuss an extension of the Standard Model with enlarged Higgs and neutrino sectors. 
Our main interest is the neutrino sector. Since we need the Higgs sector for the radiative 
neutrino masses, we give a short overview of the properties of the Higgs sector that we use 
in our calculations. 

2.1 The Higgs sector 

The authors of ref. [16] discuss the basis independent formulation of the general two-Higgs- 
doublet model (2HDM). Using their definition of the Higgs basis, we can write the two 
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complex doublets of our model in a unique way 

where the vacuum expectation value (vev) v ~ 246 GeV and the Goldstone bosons and 

appear only in the first Higgs doublet . The relations between the basis independent 
parameters defining the Higgs potential and the parameters describing the physical states 
are linear and can be easily inverted. This feature allows us to use the vev, the masses of 
the physical Higgs bosons, m^jo, rrij^o, rriffO, and mH+, and their mixing angles ai 2 and 
ai 3 as input parameters. 

The mass eigenstate for the charged Higgs boson corresponds directly to the field 
with the mass mff +, but the mass eigenstates for the neutral Higgs bosons with the masses 
mjjo, rufjo, and rrif^o, respectively, are linear superpositions of the neutral fields T-L^r-, 
and 7^21 • Following the formulation of Grimus and Lavoura [3, 4] these linear superpositions 
are conveniently expressed by 

riH ^ riH 

hi = 4 >l = V 2 Re( 6 l/) = v/ 2 ^Re( 6 ^^. 0 °) = ^ E ’ ( 2 - 2 ) 

i=i i=i 

where cjP are the neutral parts of the Higgs doublets without the vev: cj^i = (pl — vj^fi 
and 02 = 02- The ”b-vectors” are ‘Inn unit vectors bk E of dimensions tih x 1- We 
discuss those vectors in the general case in appendix A. There we also show how to obtain 
the following parametric values for the vectors b: 

ftco=f , h=( y b2=( y 63 =^ ^2.3) 

\o J y-S 12 - *ci2Si3 y y C 12 - ^sl2Sl3y y^cisy 

where cij = cosaij and sij = sinaij are given by the angles that describe the mixing of 
the neutral Higgs fields. 

Restricting ourselves to CP conserving cases we use the analysis of ref. [14] , where the 
authors discuss the CP-invariant Higgs potential in the 2HDM framework under various 
basis-independent conditions. The possible overall phase, that can be written in front of 
the second Higgs doublet and that acts like a mixing angle 023 between ?^ 2 r ^ 2 i) 

used to define the CP-property of the mass eigenstates, corresponding to their coupling to 
fermions, taking H 2 to be CP-even and to be CP-odd. This justifies the disctinction 
between the fields hi and H^. 

Having a fixed SM Higgs mass m^jo and assuming it to be smaller than the other 
two, non-degenerate neutral Higgs boson masses, we have four conditions (case I, case H, 
case Hla with rrif^o < ruj^o, and case HIb with m^jo > rrijjo), which are listed in Table 1. 
A study reported in [16] suggests that generality is not lost assuming —| ^ ai 2 ,«i 3 < f • 
We performed the numerical analysis of the neutrino mass spectrum considering the named 
cases. In some situations we refer to those cases as “scenarios.” 
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Table 1. Basis-independent conditions for a CP-conserving 2HDM scalar potential and vac¬ 
uum [14]. Uij label the mixing angles of neutral Higgses, muo and rrij^o denote the masses for 
CP-even and CP-odd Higgses. 


2.2 The Yukawa couplings 

Using the vector-and-matrix notation, the Yukawa Lagrangian for the leptons is expressed 
by [3, 4] 

nH=2 \ / \ 

KM + H-c., (2.4) 

where (f)k = The quantities Ir and ur are the vectors of the right-handed charged 

leptons and the right-handed projection of the neutrino singlets, respectively, (.r and vr 
form the lepton doublet under the weak interactions and combine with the Higgs doublets 
(j)k to form S?7(2)weak-invariant terms. They are also vectors in generation-space, ur = 
3 denoting the three generations of the SM. The Yukawa coupling matrices T^ have a 
dimension urX ur, while have a dimension ur x ur, where ur is the number of singlet 
neutrino fields. 

Taking the bilinear terms of eq. (2.4), which means taking only the vev from the Higgs 
doublets, we get the Dirac mass terms for charged leptons and neutrinos: 




= diag(me,m^,mr) 


(2.5) 


and 

Md = 

These have to be diagonalized using a singular-value decomposition (SVD) like in the 
SM to get the correct definition for the mass eigenstates that will describe the physical 
particles. Having done this transformation to the mass eigenstates, which we write down 
as the fields appearing in eq. (2.4), the respective transformation matrices reappear in two 
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unique combinations, Vckm and VrmnS) in the interactions with the charged gauge bosons 
or the charged scalar bosons and G^, giving the charged current Lagrangian 


£ 


CC 


g 

^/2 


(lI^PlVpmnsi^l + H.c. , 


(2.7) 


where g is the SU{2) gauge coupling constant. We give this part of the Lagrangian only 
as a reference, to show what neutrino experiments measure, as this PMNS matrix Lrmns 
is the basis for the interpretation of experimental data in the “3 x 3” neutrino mixing 
model [1]. 


2.3 Neutrinos at tree level 

The singlet neutrinos, added to the SM, are neutral with respect to all gauge groups of 
the SM. This offers the possibility that they are Majorana particles, allowing to write a 
Majorana mass term for them. Since the Lagrangian has to be a scalar with respect to 
Lorentz transformations, we have to combine a spinor with itself in a Lorentz invariant 
way. The only chance for Dirac spinors is to use the charge conjugation matrix C, which 
also appears in the definition of the Lorentz covariant conjugation^ for spinors 

i) := . ( 2 . 8 ) 


The Majorana condition can now be written as 

V'M = , (2.9) 

where is the Majorana phase. Assuming to be nn Majorana fermions we can write 
down a Majorana mass term as 

■^Majorana—mass — 2 t^rMrvr + H.C. = I^rMrCur + H.C. , (2-10) 

where the order of Mr and C is irrelevant, as these matrices act on different indices of the 
spinor rr: C is a 4 x 4-Dirac matrix, connecting the spinor indices of vr, whereas Mr is a 
symmetric ur x ur matrix, acting on the ’’generation” index of vr. Since the mechanism 
generating the Majorana mass is not known, we assume the singlets vr to be already in 
the mass eigenstate of Mr. This means, we assume Mr to be diagonal, containing the 
Majorana masses of the heavy singlets: Mr = Mr. 

Together with the Dirac mass, coming from the Yukawa terms eq. (2.4), the mass terms 
for the neutrinos are 


Cu- 


z/—mass 


= -vrMrI'l - \i>rMri>r + H.c. 

= —\^rMdvl — ImMJ/Ur + ^urMrCu]^ + H.c. 


I'R 


0 Ml 
Mr Mr 



+ H.C. 


( 2 . 11 ) 


^For a very clear and exhaustive description of the difference between Majorana and Dirac spinors, see 
ref. [17]. 
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and can be written in a compact form by introducing an {nL + nR) x (nL + tir) symmetric 
neutrino mass matrix 




( 0 Ml \ 

\Mn Mr) 


( 2 . 12 ) 


The neutrino mixing matrix M^, can be diagonalized [2, 4] using the properties of the 
singular-value decomposition of a symmetric matrix 


U = m = diag , (2.13) 

where m;. and m/j. are real and non-negative with mass-ordering mi^ < mi^ < mi^ and 
< • • • < order to implement the seesaw mechanism [18, 19] we assume that 

the elements of Mr are of order mu and those of Mr are of order rriR, with <C ruR. 
Then, the neutrino masses nii. with i = 1, ... ,nL (where tir = 3), are of order ml/mR, 
while the masses m/j. with i = l,n/j (where ur = 1 or 2), are of order niR. It is useful to 
decompose the (n^ -|- hr) x (ni + ur) unitary matrix U as [2, 4] 


U = 



(2.14) 


where the submatrix Ur is riRX {ur + nR) and the submatrix Ur is ur x {ur + nR). These 
submatrices obey certain unitarity relations: 


UrUI = K^, URUl = ln^, URUl = 0n^y,n^ and uIUr + UIUI 

Combining with eq. (2.13), we can obtain the following relations: 

UrTuU]^ = 0, UrttiUI^ = Mr, and URihUl = Mr . 


(2.15) 

(2.16) 


With these submatrices of U, the left- and right-handed neutrinos can be written as 
linear superpositions of the ur + ur physical Majorana neutrino fields Xi- 


UR = UrPrx and Or = U^Prx or ixr = UrPrx , (2.17) 


where Pr and Pr are the projectors of chirality. 

Switching to the physical Majorana mass states X) we have to express the couplings 
using the matrices Ur and Ur. The loop corrections, described in the next subsection, 
depend on the neutrino couplings to the Z-boson and to the neutral Higgses. Interaction 
with the Z boson is given by 


rU) 

^nc 





Pl {uIUr) 


Pr{uIui) 


X, 


(2.18) 


where is the cosine of the Weinberg angle. The Yukawa couplings for the neutral scalars 
take the form 



=4 


^ E hlx[{u],A,^UR + U]:AlU*R)PR 


2y/2 

+ (ul^Iur + UrKpI}Pr 


X- 


(2.19) 
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where we treat the Goldstone boson as /i^. The Yukawa coupling is the result 
of rewriting the Yukawa Lagrangian eq. (2.4) using the physical Higgs fields defined in 
eq. ( 2 . 2 ): 

riH 

= ^{bk)jA,. ( 2 . 20 ) 


2.4 Loop corrections to the neutrino masses 


We are interested in radiatively generated neutrino masses at one-loop level. The largest 
influence from the corrections to the neutrino mass matrix has the neutrino Majorana mass 
term 5Mi, since this submatrix is zero at tree level, = 0. The contributions from 

charge-changing currents are subdominant [3, 4]. 

We calculate the radiative light neutrino masses following ref. [3]. Once the one-loop 
corrections are taken into account the neutral fermion mass matrix is given by 




f 6 Ml Ml + 5Ml \ 
[md + 5Md Mr + 6Mr J 


f 6 Ml Ml \ 

\Md Mr) ^ 


( 2 . 21 ) 


where the 63 x 3 matrix appearing at tree level ( 2 . 12 ) is replaced by a symmetric matrix 
6 Ml- This correction dominates among all the sub-matrices of corrections. The one-loop 
corrections to 6 Ml originate via the self-energy function S'^^'^^(O) (where X = Z,G^,H^) 
that arises from the self-energy Feynman diagrams. The contributions are evaluated 

at zero external momentum squared (p^ = 0). The neutrino couplings to the Z, Higgs 
and Goldstone G^ bosons are determined by eqs. (2.18) and (2.19). Each diagram contains a 
divergent piece but the sum of the three contributions yields a finite result. The expression 
for one-loop corrections is given by (see e.g. [3]) 


5Ml = 


1 


^ 327 r 2 ^b,URm 
k=l 


m 


— 1 I In 


m 


m 




m 






m 


m 


t/^A, 




( 2 . 22 ) 


where the sum index k runs over all neutral physical Higgses 


3 Results for the case ur = 1 

3.1 General parameterization 

First we consider a minimal extension of the standard model by adding only one right- 
handed neutrino field vr to three left-handed fields vr- This simple model is useful because 
it has a small number of parameters. It allows us to obtain relations between the free 
parameters and the light neutrino masses. 

For the study using the general parameterization (reported in [7]), we fix the Higgs 
sector by choosing a specific GP-conserving set of vectors h: 
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We use a parameterization of the Yukawa matrices Ai and A 2 in the following form: 


A. 


%/2 


■mo 

- a, , 


i = 1,2. 


(3.2) 


We further assume that the linearly independent vectors are normalized, |ai| = |a 2 | = 1; 
and equal total strength of the couplings. Using the block form, diagonalization of the 
symmetric neutrino mass matrix at tree level (2-12) can be written as 


= uL. 


f 03x3 rnnai 
y mnaj Mr 


Utrpo — 


V 0 WV' 


(3.3) 


The non-zero masses in and can be determined analytically by finding the eigen¬ 
values of the hermitian matrix These eigenvalues are squares of the neutrino 

masses, = diag(0, 0, and The solutions 


2 (0) (0) 

rriR =mj^ m] , 

Y _ / (0) _(0) 


m 


r )‘ 


(3.4) 

(3.5) 


correspond to the seesaw mechanism. 

We diagonalize the tree-level neutrino mass matrix using a diagonalization matrix 
Utree made of two diagonal matrices of phases and three rotation matrices: 


Utree = U^(<(-i)Ui2(ai)C/23(a2)U34(/3)Ui, (3.6) 

The angle j3 is determined by the masses and tan^(/3) = The 

dependency of 4>i and ai on rriD and oi may be expressed analytically. The diagonalization 
matrix Utree is unitary because the rotation matrices Uij and the diagonal phase matrices 
and Ui are all unitary. 

Diagonalization of the nentrino mass matrix with the one-loop corrections included, 
in eq. (2.21), is performed numerically using a unitary matrix 

Uloop ~ UegvU(p((/7l, V^2) V^s)) 

where Uggv is an eigenmatrix of mI)^\ and is a phase absorption matrix. As 
discussed already in ref. [2], the heaviest light neutrino obtains mass at tree level from 
the seesaw mechanism. The second light neutrino obtains mass from radiative corrections. 
The lightest neutrino remains massless. 

The numerical evaluation of the model parameters and of the masses of the light 
neutrinos proceeds in several steps. First, the mass matrix for the tree level is constructed. 
The lightest neutrino remains massless in this model, mi^ = 0. Using the central values 
of the experimental nentrino mass differences, we estimate the mass of the heaviest light 
neutrino and use it as an input parameter m™. The value of niR is another input parameter. 
Using the seesaw relations (3.4) and (3.5), we evaluate mo = mRim}^. The vector oi 
is generated randomly. This fully determines the tree-level nentrino mass matrix M^\ 
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Solving the eigenvalue equation (3.3) we obtain the tree-level neutrino masses mf‘\ 
and the diagonalization matrix 17tree- The next step is to evaluate the one-loop corrections 
to the neutrino mass matrix 6 Ml. Additional input parameters of 0,2 enter the procedure. 
Diagonalization of the corrected neutrino mass matrix yields masses for two light 

neutrinos. If the obtained mass differences are compatible with the experimental data on 
neutrino oscillations, the model parameter set is kept. Otherwise, another set of input 
parameters {m/?, oi, 02 , } is generated. 

The study suggested a lower limit of 830 GeV for the mass of the heavy singlet [7]. 

3.2 Reduced parameterization 

For the Higgs sector we use the values of the orthogonal complex vectors h listed in Table 1. 
The mass of the lightest neutral Higgs is fixed at = 125 GeV. The masses of heavier 
neutral Higgses rrif^o and are generated randomly in the range from 126 to 3000 GeV. 

The light neutrino helds can be transformed in such a way that = (0,0, l), and 
aj = ( 0 , n, n') with real numbers mD,n > 0, and a complex number n'. Due to the assumed 
normalisation conditions |ai| = |a 2 | = Ij there are only two independent parameters, 
namely, the real number n (n < 1), and a complex phase </>: 

37 =(0,0,1), (3.8) 

aj =(0, n, 1 — n^) . (3-9) 

A similar case was studied in sect. 4 of ref. [2] without the assumption |a 2 | = 1. Since the 
goal was to demonstrate that one of the massless neutrinos can obtain mass through the 
1-loop radiative corrections, only the expressions are given there. In addition, the masses 
of the lightest and the heaviest light neutrinos do not change going from the tree level to 
the 1-loop-corrected level in the analysis of ref. [2]. 

The numerical evaluation of the model parameters and the masses of the light neutrinos 
in the case of the reduced parameterization is done similarly to the general case, described 
above. At the one-loop accuracy, our model will predict a vanishing mass for the lightest 
neutrino. If we use this value directly, the measured neutrino mass differences would lead to 
highly restricted values of the neutrino masses. In order to reduce the impact of a starting 
point to our analysis we allow the lightest neutrino to have a small non-vanishing mass. 
Using the central values of the experimental neutrino mass differences, we determine 
the largest initial value of the light neutrino masses mj“. Selecting the value of mu, the 
value of mo is evaluated from the seesaw relation, eq. (3.4). Since the vector oi is fixed in 
this case, the tree-level neutrino mass matrix is fully defined as an input quantity at 
this point. 

Solving the eigenvalue equation, we get the tree-level masses of the heaviest light neu¬ 
trino and of the heavy neutrino . The diagonalization matrix I/tree is constructed 
from a rotation matrix and a diagonal matrix of phases: Utree = Then one-loop 

corrections to the neutrino mass matrix are evaluated, and the parameters defining 02 enter 
into the further evaluation. Diagonalization of the 1-loop neutrino mass matrix is 

performed numerically with a unitary matrix Uioop as in the general parameterization case. 
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Figure 1. (Color online) The 3D histograms of 023 = 0atm oscillation angles for ur = 1. The left 
plot represents the case II of the Table 1. The right plot shows the case I. The filled boxes indicate 
the 3 ct experimental boundaries [12], the blue vertical lines denote the experimental central value 
of 023- 

After this procedure two light neutrinos have masses and For the calculation 

of ^rnfj we assume the mass of the lightest neutrino to be the input value with the 
justification, that it could be generated by a two-loop contribution. If the calculated mass 
differences 

^("^ 12 ^ = and A(m^ 3 y = (3.10) 

and the determined value of 023 ~ 0atm (as described in the next paragraph) are compatible 
with the experimental data on oscillations, the model parameter set is kept. Otherwise, 
another set of input parameters {m^”, n, (j), aij, rrijjo } is generated. 

The results of the 1-loop corrected calculations are subject to the constraints from 
the experimental data on solar and atmospheric neutrino oscillations [12]. The neutrino 
mixing angles are determined from a factorization of the diagonalization matrix 17ioop into 
terms where the PMNS matrix is included explicitly, following the ideas of ref. [13]. (The 
method is described in detail in appendix B.) For the case ur = 1 it is possible to find exact 
analytical expressions for the mixing angles, Dirac and Majorana as well as non-physical 
phases. The oscillation angles are constrained by the experimental data in the sense that 
having a randomly generated set of the input parameters we derive the 1-loop corrected 
results and require that the estimated mixing angles are consistent with the experimental 
values in the 3 cj range. It should be noted that this reduced parameterisation has only one 
non-vanishing neutrino mixing angle, namely, the atmospheric neutrino mixing angle 0atm, 
because the vectors ai and 3,2 have the same vanishing component (an = 021 =0). In a 
more general case, all three oscillation angles are non-zero. 

Distributions of the atmospheric oscillation angles for the cases I and II are shown in 
figure 1. The plot on the right shows the case I that is similar to the cases Ilia and Illb. 
This distribution is narrower than the one obtained in the case II. The distribution for the 
case I has a well pronounced peak at around 45 degrees, while the probability to have an 
oscillation angle 023 in the range from 25 to 65 degrees is nearly flat in the case II. 
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Figure 2. (Color online) Calculated masses of two light neutrinos as a function of the heaviest 
neutrino mass mh for 3 different scenarios of Table 1 and nu = 1. The black solid line represents 
the initial parameter m-in- Values of mi^ and mi^ are shown as a band where the color corresponds 
to the mass of the heavy neutrino mh, given on the horizontal axis. This color code is also used in 
hgures 3 and 4. 


Since only two light neutrinos acquire mass in the case of ur = 1 with 1-loop correction 
included, only normal ordering of neutrino masses is possible. Assuming that the lightest 
neutrino mass mq = 0, we obtain a hxed spectrum of light neutrinos with mq = 8.7 ± 
0.3 meV and mq = 50.6 ± 2 meV. However, the peak of the derived distribution of the 
oscillation angle 0atm is shifted from the experimental value. 

Figure 2 illustrates the neutrino mass spectrum for three different scenarios given in 
Table 1, assuming ttir = m\^ ^ 0. All scenarios with hr = 1 produce distributions of 
neutrino masses which differ from each other by value. In the hr = 2 case (which is 
discussed below), all distributions look very similar, regardless of the scenario. The cases 
Ilia and Illb have similar distributions, therefore no reference to the case ’a’ or ’b’ is given 
in the plot. The mass of the heaviest neutrino mq reaches the highest value of 140 meV, 
when rrih = 10^ GeV in the case III. The lowest value of mi^ = 70 meV is obtained when 
rrih = 10^^ GeV in the cases I and II. The mass of the intermediate light neutrino mq 
reaches the highest value of 121 meV, when rrih = 10“^ GeV in the case III, and the lowest 
value of mq = 48 meV is reached when ruh = 10^*^ GeV in the case I. 

We observe a particular relationship between the values of the masses of light and heavy 
neutrinos. The calculated masses of the light neutrinos decrease, if the heavy neutrino mass 
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Figure 3. (Color online) Values of the neutral Higgs masses 171^0 and m^o as a function of the 
heaviest neutrino mass ruh for 4 different cases of Table 1 in the n/j = 1 model. The scale of the 
ruh values is shown on the right. The mass of the SM Higgs boson is fixed to rrijjo = 125 GeV. 


increases. This dependence emerges from the relation of ruh to 0atm- The value of min has 
to get lower as vr gets heavier in order to keep the most probable value of the oscillation 
angle within the experimental range of 3(7 (this range is marked by boxes in figure 1). The 
narrow bands of the values of mi^ and m/g are formed by a relatively high value of min and 
the restrictions on the mass differences. The results suggest that the lowest limit of the 
heavy neutrino mass is 10^ GeV. 

The allowed values of the Higgs masses (other than the SM Higgs) are illustrated in 
figure 3 as a function of the heavy singlet mass m/j. A band structure is formed according 
to the choice of vectors h (see Table 1) and the values of the free parameters n and which 
are displayed in figure 4. When the mass of the heavy singlet rrih is increasing, the values 
of the Higgs masses tend to decrease in the 2HDM model cases I and H. This tendency 
is absent in the case HI, where the allowed values of form a narrow band. Larger 

values of rrih make the masses of more similar, although their difference does not 

disappear. The range of the allowed masses of the heavier Higgs boson starts at 500 GeV. 

The values of the Higgs masses displayed in figure 3 satisfy the experimental restric¬ 
tions [1] of the oblique parameters S, T, and U, introduced by Peskin and Takeuchi [20]. 
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Figure 4. (Color online) Values of the free parameter 023 as a function of the heaviest neutrino 
mass mh for the cases of Table 1 and nn = 1. The scale of the rrih values is shown on the right. 


These oblique parameters define combinations of observables that quantify deviations from 
the SM predictions. The experimental electroweak precision data yields values compatible 
with SM. We estimated these parameters using the algorithm implemented in the two Higgs 
doublet model calculator (2HDMC) [21] taking recent values of SM parameters [1]. The 
values of the oblique parameters are functions of masses of the three neutral (rrifjo, rn^jo, 
and rrifjo) and one charged {mH±) Higgs boson, and an additional factor that is related 
to the mixing angles of the neutral Higgses. This factor, defined as sin(/3 — a) in ref. [21], 
equals to cos(ai 3 ), cos(q;i 2 ), and 1 in the scenarios I, H, and HI of Table 1, respectively. 
The estimated masses of the neutral Higgses and their mixing angles in our model are 
discussed above. The mass of the charged Higgs mH± is allowed to take values in the range 
of 126-3000 GeV. This defines the parameter space to estimate the values of S, T, and U. 
The calculated curves of S and U are mostly contained within the experimental bounds. 
The experimental limits on T are narrower as compared to the calculated distribution of 
T. A similar observation on the restrictive power of the oblique parameters is made in 
ref. [22]. These limits can be used to restrict the mass of the charged Higgs boson. If the 
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values of all three oblique parameters are compatible with the experimental data and the 
allowed value of mH± is within the studied range, the tested set of model parameters is 
kept. Otherwise, another set of parameters is generated. Using a stronger constraint of 
U = 0, which is also compatible with SM, about 35%-45% of the values in figure 3 satisfy 
the experimental restrictions [1] . In this case the majority of rrifjo values fall in the range 
of 126-1500 GeV. 

Distributions of the free parameter 023 , part of the neutrino couplings to the second 
Higgs doublet, are shown in Fig. 4. When the mass of the heavy singlet nih increases, 
the absolute value of the parameter |a 23 | also increases. In case I the real part of 023 
varies between —1 and 1 and the imaginary part varies from —0.18 to 0.18. In case II 
real and imaginary parts behave like interchanged: the real part of 023 varies between 
—0.18 and 0.18 while imaginary part varies from —1 to 1. In case Ill-a the real part of 
023 varies between —0.34 and 0.34 while the imaginary part is restricted to the intervals 
(0.04,0.17) and (-0.17,-0.04). In case Ill-b real and imaginary parts are exchanged like 
between cases I and II: the real part of 023 is restricted to the intervals (0.04,0.17) and 
(—0.17, —0.04) while the imaginary part varies between —0.34 and 0.34. 

In summary, tuning the value of min, and restricting the light neutrino mass differences 
to the experimental central values within Icr and the angle of oscillations 623 = Oatm within 
3a, we determined the lowest limit of 10^ GeV for the mass of the heavy neutrino singlet. 


4 Results for the case ur = 2 


4.1 General parameterization 

If we add two singlet fields vr to three left-handed neutrino fields i'l, the radiative correc¬ 
tions give masses to all three light neutrinos. In the general case we parameterize 


A* 


^ f moj] \ 
V y moSl ) 


(4.1) 


with 12 complex parameters of Oj and bi, where i = 1,2. We further assume that the 
vectors are normalized, |aj| = | 6 j| = 1 . A specific set of the vector b values (3.1) was used 
for this study. 

Numerical evaluation of the model parameters and the masses of the light neutrinos 
is performed in several steps. First, the neutrino mass matrix for tree level is constructed. 
The lightest neutrino is massless at tree level, = 0. Taking = 0, the masses of 
the other two light neutrinos, and mj“, are estimated from the experimental neutrino 
mass differences like for the case ur = 1. Entries of the heavy neutrino mass matrix, 
Mr = diag(mRi , m-Rj), are input parameters. The eigenvalue equation in the block form 
is: 


= uL 


^tree 


/ 03x3 rnDaai mojii 

\ 


M. 


Ui 


tree 


R 


(0 ^ 
I 0 A/f j 


(4.2) 
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2HDM 


scenario 


Normal hierarchy 


A 


1 


Inverted hierarchy 


A 


1 


011 012 013 
0 6i2 bi 3 


0 022 023 

0 b 22 623 


Oil O12 Oi 3 

0 612 bi 3 


021 022 0 
b 2 i 0 623 


II 


011 012 013 
0 612 bi 3 


021 022 0 
0 622 b23 


011 012 013 

0 612 bi 3 


0 022 023 
&21 0 623 


IIia,b 


Oil 012 0 
0 612 bi 3 


0 O 22 023 
&21 0 623 


Oil O12 Oi 3 

0 612 bis 


02 I 022 0 
b 2 i 0 623 


Table 2. Textures of Dirac matrices used for calculations of the light neutrino mass spectra for 
normal and inverted hierarchies in four different cases of the 2HDM vectors b, listed in Table 1. 


where = diag(m|||\ with The values ofrriDa ond rriD^ are evaluated 

through the seesaw equations, relating niD, mn on one side and mi, on the other side: 

^Da = (4-3) 

mo, = 

m\. ~ * = 1,2. (4.5) 

The diagonalization matrix for tree level C/tree = is composed of the eigen- 

matrix of (denoted by and a diagonal phase matrix C/^. At the second 

step we evaluate one-loop corrections to the neutrino mass matrix. Diagonalization of 
is performed with a unitary matrix Dioop = Ulgv^U,p{(pi), where C/igv’’ is the eigenmatrix 
of and 17^, is a phase matrix. This procedure yields masses for all three light 

neutrinos. If the calculated mass differences are compatible with the experimental data 
on neutrino oscillations, the model parameter set kept. Otherwise, another set of input 
parameters {mji^ 2 J^ijbi,mjjo is generated. 

The numerical analysis with the full set of parameters, constrained only to the exper¬ 
imental mass differences of the light neutrinos, suggests that heavy singlets should have 
mass greater than 100 GeV [7]. 


4.2 Reduced parameterization 

Studying the influence of the randomly generated parameters we found that a reduced 
number of parameters is sufficient to fulfil the experimental criteria [12] of Atuq, 

012 , 013 and 023- We selected several “textures” of Dirac matrices (i.e. the patterns of 
non-zero components), which allow the most accurate agreement to the experimental data. 
The texture of the matrix Ai has the largest impact to the results of oscillation angles. 
There are 3 best versions of the Ai textures: 


an ai 2 013 
0 bi 2 bi 3 


I ail 0 ai3 

y 0 612 613 


and 


( 0,11 ai 2 0 
y 0 612 bi 3 


(4.6) 
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The textures of the matrix A 2 play a subdominant role for the results. In our calculations 
we tailored the textures to the 2HDM scenarios. As a result, the second texture of Ai, 
having 012 = 611 = 0, was not used. The applied textures of Ai and A 2 are listed in 
Table 2. All non-vanishing components can have complex values. 

Studies of textures with one or two zero entries in the models with two heavy singlets 
have been reported in refs. [23-28]. Neither of those studies considered the SM exten¬ 
sion by 2HDM. We obtain experimentally compatible neutrino oscillation values as in the 
mentioned references, but our textures are tuned to provide the model parameters in best 
agreement with the experimental data. It is worth mentioning that the textures listed 
in Table 2 may be a natural consequence of a certain flavor symmetry, as discussed in 
refs. [29-31]. 

The numerical evaluation of the model parameters and of the light neutrino masses, 
using the textures listed in Table 2, is done similarly to the general case. Having picked 
a mass for the lightest neutrino as an input value m\n = , the masses of the other two 

light neutrinos, and r?T,]“, are estimated from the central values of experimental neu¬ 
trino mass differences. Entries of the heavy neutrino mass matrix, = diag(mK^, m^^), 
are input parameters. The seesaw relations (4.3) give the values of moa mof,- The 
tree-level neutrino mass matrix is diagonalized, and Utvee is obtained. Then one-loop 
corrections to the neutrino mass matrix are calculated, and the parameters defining A 2 
enter into the further evaluation. Diagonalization of the corrected neutrino mass matrix 
yields masses for all three light neutrinos. If the calculated mass differences are 
compatible with the experimental data on oscillations, including the three neutrino mix¬ 
ing angles, the model parameter set is kept. Otherwise, another set of input parameters 
di,bi, aij,mjjo J is generated. 

To increase the efficiency of random sampling, we chose the initial mass rriin in the 
range of 16-20 meV according to the agreement of the calculated results to the experimental 
values of the oscillation angles. Unlike the case ur = 1, the oscillation angles are less 
affected by the value of m[^. Calculation of radiative corrections is done for the scenarios 
of the orthogonal complex vectors b listed in Table 1. The mass of the lightest neutral 
Higgs is fixed at Tn^jo = 125 GeV but the masses of heavier neutral Higgses, and 
rrif^o, are generated randomly in the range from 126 to 3000 GeV. 

Figure 5 shows the distributions of the oscillation angles for the scenario HI using 
the best textures for Ai and A 2 . Both normal and inverted neutrino mass hierarchies 
are shown. In the case of the inverted hierarchy the peaks of the distributions agree very 
well with the experimental bounds of all oscillation angles. In the case of the normal 
hierarchy, the most probable value of 612 is slightly different from the experimental value. 
It is around 55 degrees, instead of the expected 33.4 degrees. However, if we extract 612 
from sin^( 20 i 2 ) we will hnd two solutions in the first quadrant. Following PDG [1], the 
experimental boundaries are in the regions 32.6-34.3 degrees and 55.7-57.4 degrees which 
agree with the peak of the atmospheric angle. In hgure 5 the resulting values of ^13 are 
more localized as compared to the inverted hierarchy case. The distributions of 623 agree 
with the experimental values in both hierarchies. The same tendency of the mixing angle 
distributions is seen in all three scenarios of the 2HDM basis. It should be noted that for 
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Figure 5. (Color online) The 3D histograms of the oscillation angles for the scenario III in the 
case nn = 2. The plots on the left (right) correspond to the normal (inverted) hierarchy of the light 
neutrino masses. Filled boxes mark the applied experimental boundaries of 3 ct, blue vertical lines 
mark the experimental central values. 


the general parameterization case (where we have 12 complex parameters) all peaks of the 
oscillation angle distributions take values at approximately 45 degrees. 

The neutrino mass spectrum is analyzed assuming that the masses of the heavy neutri¬ 
nos are nearly equal, ruh^ ^ (the ratio is 0.999). This simplifies the analysis 

and does not change the distributions of the light neutrino mass spectra in general terms 
(we can compare figure 6 to figure 2 in ref. [7]). The masses of the light neutrinos mi^ 
are estimated when the masses j vary in the range from 500 GeV to 10^® GeV, see 
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Figure 6. (Color online) Masses mi. of the light neutrinos as functions of the heaviest neutrino 
mass, max(m;n, 771/12), for the scenarios of 2HDM in the case of tir = 2 (the scenarios Ilia and 
mb are very similar). The plots on the left represent the normal hierarchy, the plots on the right 
represent the inverted hierarchy of the light neutrinos. Wide bands indicate the area of the most 
frequent values of the scatter data. The nearly-degenerate masses mi^ and mi^ are shown separately 
in the lower right plots for the inverted hierarchy. 


figure 6. Both normal and inverted hierarchies are shown. The dependency of the light 
neutrino masses on the values of the heavy neutrino masses is similar in all three scenar¬ 
ios. If normal hierarchy is assumed, the lightest neutrino mass that is generated by 
the one-loop corrections varies from O.OI to 20 meV. The most frequent values lie around 
3 meV, when the mass of the heavy neutrinos is ruh^ ^ > 3000 GeV. The mass mi^ varies 
from 8 to 25 meV with the most frequent values at 10 meV. The mass of the heaviest 
light neutrino m/g is around 50 meV. If the inverted hierarchy is assumed, the range of 
values is wider and varies from 0.001 to 40 meV. The most frequent values increase and 
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Figure 7. (Color online) Values of the free parameters rrij^o and rrifjo as functions of the heaviest 
neutrino mass max(m/jj, for the scenarios of 2HDM in the case ur = 2. The plots on the left 
represent the normal hierarchy, the plots on the right represent the inverted hierarchy of the light 
neutrinos. The scale of the max{mki,rnh 2 ) values is shown on the right. 


take values of 0.5-4 meV, depending on the masses of the heavy neutrinos. The values of 
m /2 and m/g are nearly degenerate and vary from 47 to 55 meV. The most frequent values 
are in the range of 48-51 meV. 

Figure 7 illustrates the allowed values of the Higgs masses depending on the masses of 
the heavy singlets m/^^ g • Scenarios I and II are rather similar in dependencies, namely, an 
increase of the heavy singlet mass leads to the decrease of the Higgs masses. The mass of 
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the second Higgs boson tends to be different from the mass of the third Higgs boson. The 
scenario HI has different dependencies. The heavy Higgs masses tend to be equal for large 


values of the heavy singlet masses, and tend to be independent of it. 

The values of the Higgs masses displayed in figure 7 also satisfy the experimental 
restrictions [1] of the oblique parameters S, T, and U, as discussed in previous chapter. 
Using the stronger constraint of C/ = 0, about 40%-65% of the values shown in figure 7 
satisfy the experimental restrictions. In this case the majority of m^o values fall in the 
range of 126-1000 GeV. 

The presented results are obtained using the tuned textures of Table 2. An alternative 
method to reduce the parameter space of vectors a and b could be used. For example, 
instead of setting an entire component {aij or bij) to zero, the parameter space could be 
limited restricting the values of these vectors to real numbers. A study of these textures 
will be reported in the future. 

The most reduced parameterization of Dirac matrices, that still allows experimentally- 
compatible results of the light neutrino mass differences, has only 4 independent real pa¬ 
rameters: 



(4.7) 


(4.8) 


with |aj| < 1 and \bi\ < 1, i = 1,2. The neutrino oscillation angles evaluated in this 
strongly-reduced parameterization do not have the most-probable values in the experimen¬ 
tally determined range. 

5 Summary 

The seesaw mechanism is one of the most successful extensions of the SM which explains 
neutrino masses. Finite corrections to the neutrino mass matrix arise from one-loop dia¬ 
grams mediated by a heavy neutrino. In our model the Higgs sector is constructed with 
two Higgs doublets and a CP-invariant Higgs potential which allows to distinguish four 
conditions for vectors b and thereby determine the scenarios (see Table 1) for numerical 
calculations. The SM Higgs mass is fixed to 125 GeV. By diagonalizing the neutrino mass 
matrix we obtain light neutrino masses and derive their oscillation angles. Sets of free 
parameters have been selected according to the distributions of oscillation angles within 
the experimental boundaries, when also the masses of the light neutrinos give the measured 
neutrino mass differences. In this paper we have studied two cases when one or two heavy 
neutrinos are added to the three light neutrinos. We refer to those cases as nn = 1 and 
Ur = 2 . 

In the Ur = 1 case with a general parametrization of the Dirac matrices there are six 
free complex parameters. The numerical analysis suggests a lower limit of 830 GeV for the 
heavy singlet mass. However, the large number of free parameters makes it difficult to find 
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the correlations among them. By reducing the number of parameters we can study relations 
between them and the dependency on the heavy neutrino mass. The minimal reduction of 
free parameters, namely, = ( 0 , 012 , 013 ) and oj = ( 021 , 022 , 023 ), allows to estimate all 
three oscillation angles 0 i 2 , ^13 and 023- The peaks of their distributions agree well with 
the experimental bounds. We presented results of a strongly reduced parameterization, 
aj = ( 0 , 0 , 1 ) and oj = ( 0 ,n,e*'^\/T — n^). According to the chosen minimal set of the 
free parameters only the angle 623 can be estimated. The calculated masses of the light 
neutrinos decrease, when the heavy neutrino mass increases. This dependency emerges 
from the relation of mh and 023- Tuning the initial value min, and restricting the light 
neutrino mass differences to the experimental central values and the oscillation angle 623 
within 3cj, we determined for the mass of the heavy neutrino singlet a lowest limit of 
10^ GeV. When the mass of the heavy singlet nih is increasing in the scenarios I and II of 
the 2HDM model, the allowed values of the non SM Higgs masses tend to decrease. This 
tendency is absent in the scenario III, where the Higgs boson masses m^o get closer to 
each other, but stay different. We find a lower limit for the allowed values of the heavier 
Higgs boson mass of about 500 GeV. The values of the free parameters depend weakly on 
the mass of the heavy singlet mh- 

The general parametrization of the Dirac matrices in the ur = 2 case has twelve 
complex parameters. The numerical analysis shows that the heavy singlets should have 
masses greater than 100 GeV. However, the most probable values of the neutrino oscillation 
angles are not in the experimentally determined range. We selected several textures of Dirac 
matrices, which allow the most accurate agreement to the experimental data. The texture 
of the matrix Ai has the largest impact on the values of the oscillation angles while the 
textures of the matrix A 2 play a sub-dominant role. The used textures for normal and 
inverted neutrino mass hierarchies are listed in Table 2. The neutrino mass spectrum is 
analyzed assuming that the masses of the heavy neutrinos are nearly equal, mh^ ~ 

(the ratio mh^/mh 2 is 0.999). The masses of the light neutrinos are estimated when the 
masses mh^ 2 greater than 500 GeV. The dependency of the light neutrino masses on 
the values of the heavy neutrino masses is similar in all three scenarios. An increase of the 
heavy singlet mass leads to the decrease of the non SM Higgs masses in the scenarios I and 
H. This tendency is absent in the scenario HI, where the Higgs boson masses tend 

to be equal as the masses of the heavy singlets increase. The allowed values of mj^o can 
sample the entire range. Due to the large number of free parameters it is difficult to find 
correlations among them. 

Our analysis has shown that the radiative corrections are quite sizeable and play an 
important role. They should be taken into account in the studies of the see-saw models. 
The studied case with one heavy singlet is a ’’toy” model because it is strongly restricted 
and does not provide all physical quantities. For example, the mass of the lightest neutrino 
is equal to zero, and we can evaluate only one oscillation angle. However, using this model 
it is possible to make some generalisations about the distributions of the Higgs masses and 
the correlations between the free parameters for models with a larger number of heavy 
singlets. The hr = 2 case allows the calculation of all three masses of the light neutrinos 
with a reduced number of free parameters. The finding of textures which allow the most 
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accurate agreement of the oscillation angles to the experimental data could motivate some 
future models, for example, those based on the Abelian family symmetry or another discrete 
symmetry. 


A Neutral Higgs mass eigenfields 


Some features of formalism for the scalar sector of the multi-Higgs-doublet SM is given in 
ref. [2, 4]. Here we discuss the properties of the vectors b and give expressions for their 
calculation in the case of two Higgs doublets. 

The physical neutral scalar mass eigenfields are expressed as 

riH ^ nn 

< = ^ E > (A-i) 

j=i i=i 

which are characterized by 2nH unit vectors hk G of dimensions nn x 1. In the 

matrix-vector notation, these eigenfields can be written as = y/2Yie.{h\(fP). 

The orthonormality equations for the vectors are 

nu riH 

^ {Re{bkj)Re{bk'j) + lm{bkj)lm{bk>j)) = '^Re{bljbk'j) = Sb^b^,; 
i=i i=i 

'^Re{bkj)Re{bkj') = '^lm{bkj)lTa{bkjf) = 
k=l k=l 

277.J/ 277.J/ 

^ Re(6fcj)Im(6fcy) = ^ bkjbkf = 0. 

^=1 k=l 

The vectors b^ and b^/ indicate two different states and ^, and indices j and j' 
indicate two different components of the vectors b. 

The neutral Goldstone boson corresponds to the vector with the com- 

f 2 2 2^^^^ 

ponents {bQo)j = ivj/v [2, 4], where v = ( Inif -|- \v 2 \ + • • • + ) = ‘2mw/g- In the 

case of only two Higgs doublets, and due to the rotation of the Higgs fields to make the 
vacuum expectation value a feature of the SM Higgs field, the vector bQO equals 


(A.2) 

(A.3) 

(A.4) 


bQO = 



(A.5) 


Physical Higgs fields (p^^^QO must be orthogonal to the Goldstone field which follows 
from (A.2). This leads to the condition 

^ ) = “ I] blj) = 0. (A.6) 

i=i ^ ^ i=i i=i 

To study the unit vectors b, introduced in eq. (A.l) (which is the same as eq. (2.2) in 
the text) and corresponding to the Higgs fields other than the Goldstone boson G^, lets 
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define them in the following form: 


bi = 





(A.7) 


From the orthogonality relations (A. 2 - A. 4) and due to the fixed value of bQO (A.5) it 
is possible to write the orthogonality equations for vector components in the following 
manner: 


^ 11 , ^21, ^31 £ fR; ^ 12 , b22, bz2 £ C; 

(A. 8 ) 

bh + \bk2\^ = 1 ; 

(A.9) 

bkibk'i + Re {b*k 2 bk' 2 ) = 0; 

Q Q 

(A.IO) 

0 0 

'^^k2 = '^ bkibk2 = 0 ; 
k=l k=l 

3 3 3 

(A.ll) 

= (''<=2)1" = E ('>‘2)1" = 

k=l k=l k=l 

(A.12) 


By choosing 631 , 621 and Re ( 632 ) as input variables, it is possible to express the other 
components of the vectors b by those variables by solving the equations (A .8 - A. 12 ). 
Introducing three sign-parameters S 32 im, sn, and S 22 (they can take values ±1), we can 
write 


Im(fe32) =S 32 im 1 - - [Re (632)]^ ; 

bn =sn\/l - &31 - ^21 ; 

_63i62iRe (632) + S22 l&iil |Im (632)! 


^comb 


^31 1 


P 22 = 


-Sg(63i)Sg(62i)Sg(Im(632)), if |Re(632)| < 
S 22 Sg(Re( 632 ))Sg(Im( 632 )), otherwise ; 

b22 =bcomb + W22\l'^ “ “ &co 


1 “ ^21 


■^comb ’ 


bl2= — -r~ (^31^32 + ^ 21 ^ 22 ) • 
On 


(A.13) 

(A.14) 

(A.15) 

(A.16) 


(A.17) 
(A.18) 

We introduced two intermediate parameters 6 comb and p 22 -, and Sg(x) is the sign function 


Sg(x) 


—1, X < 0, 

1, X ^ 0. 


(A.19) 


It is worth mentioning that the solutions for the parameter values, given by the equations 
(A.13 - A.17), were obtained assuming 621,^31 / il. According to the orthogonality 
relations (A .8 - A. 12 ) the free scale parameters vary in the following ranges: I 631 I < 1, 
I 621 I < y^l — 63 ]^, and |Re( 632 )| < 1 — 633 . The extreme values of ±1 for the parameters 
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621 and 631 could be obtained by the index permutation of the vectors bk (for example, 
621 = 1 can be obtained by swapping the values of bu = 1 and 612 with those of 621 and 

b22)- 

Equations (A.13 - A.17) give 8 different solutions for the vectors b, corresponding to 
two possible values of the sign-parameters Sx {x = 32im, 11, and 22). However, due to the 
structure of the one-loop corrections (2.22), only 4 different solutions of those equations 
are important, since the sign of 611 (i.e. the value of sn) does not change the values of the 
light neutrino masses. 

The expressions of eqs. (A. 13 - A. 17) are significantly simpler, if some input parameters 
are equal to zero. This can lead to further simplification after introducing trigonometric 
functions. Let us study the case, when Re ( 632 ) = 0. Defining 631 = sin(ai 3 ), 621 = 
sin(ai 2 ) cos(ai 3 ), and taking S 32 im = sii = S 22 = 1 , we obtain the following parametric 
values of the vectors b: 

(A.20) 

\o J y-S 12 - *ci2Si3 y y C 12 - ^sl2Sl3y y^cl3y 

where Cij = cos{aij) and Sij = sin(ajj). 


B Neutrino oscillation angles 

Neutrino oscillation angles are introduced using the tree-level neutrino mass diagonalization 
matrix Hioop and factorizing it to contain the ordinary Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix. We introduce the formalism by discussing the 3x3 
neutrino mixing case, where the relationships are simpler. Then we discuss the cases that 
are used in current paper, namely, (3-1-1) x (3 -|- 1) and (3-1-2) x (3 -|- 2) neutrino mixing. 

The simplest case (3 x 3) considers only the SM neutrinos. It is discussed in ref. [32] 
in a slightly different notation of the matrix elements. Factorization of the rotation matrix 
with the PMNS matrix included explicitly in the case 3 -|- 3 is discussed in ref. [13]. Here 
we give formulas for the intermediate cases. 

The neutrino masses and the mixing angles are predicted from a given neutrino mass 
matrix (the “top-down” method). Exact analytical expressions for the mixing angles, Dirac 
and Majorana phases and formulas for the non-physical phases can be given for the 3- and 
4-dimensional cases. Only numerical solutions are possible in the case of 2 additional 
neutrinos (the 5-dimensional case). 


The 3-dimensional case 

First we parameterize the neutrino diagonalisation matrix by including the PMNS 
mixing matrix for the 3x3 mixing. The neutrino mass matrix can be diagonalised by 
a unitary transformation U, obtained by the singular value decomposition method, see 
eq. (2.13). Lets denote the matrix elements in the following way: 


jj{3x3) 


^ Xi X2 X3 

yi y2 ys 

y^ii 22 Z 3 


(B.l) 
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This matrix could be factorized into three terms 


C/(3"3) = [>(3) . UpuNS • 

where ?7pmns is the standard PMNS mixing matrix for Dirac neutrinos: 


b^PMNS = 



/ Ci2 Si2 0 

-Sl2 Ci2 0 
0 0 1 


C12C13 


C 13 S 12 


=13 


—C23S12 — C12S13S23 C12C23 — S12S13S23 C13S23 
\ S12S23 — C12C23S13 —C23S12S13 — C12S23 C13C23 


(B.2) 


(B.3) 


We used abreviations c^- = cos 9ij and s^- = sin0jj, where Oij and 5ij are the rotation 
angle and the phase angle, respectively. 

The two diagonal phase matrices are defined as 



/ 0 

0 \ 


II 

0 

0 , 

(B.4) 


^ 0 0 

J 



(^ ° 

0 \ 


= 

0 

0 

(B.5) 


^0 0 

glK2/2 / 


There are 9 parameters: 3 mixing angles ( 0 i 2 

iM 

CO 

1 Dirac phase ^ 13 ; 2 Majorana 


- (3) 

phases ki and K2] and the matrix ^ containing 3 non-physical and unmeasurable phases 
4>i {i = 1,2,3). 

Comparing eqs. (B.l) and (B.2) we can find the relations between the elements of the 
rotation matrix in a general form and in the factorized form: 


6*13 


arcsin (|x 3 |), 623 = arcsin 



012 = arcsin 



(B. 6 ) 


y = arg(x 2 ) - arg(xi), y = arg(a: 3 ) - arg(xi) + (^ 13 , (B.7) 

(/)i = arg(xi), (/)2 = arg(xi) - arg(x 3 ) + arg(i/ 3 ) - ( 5 i 3 , (B.8) 

</>3 = arg(a:i) - arg(x 3 ) -h arg(z 3 ) - 613, (B.9) 

613 = arg(x2) - arg(x3) + arg(y3) + i\ni - l^sl ) + 3 : 23 : 31 / 3 ^ ^ (B.IO) 

\ \Xl\\Z3\ J 

These relations are obtained comparing eq. (B.2) with the following matrix elements from 
eq. (B.l), forming the upper-triangular matrix: xi, X 2 , X 3 , 1 / 2 , 2 / 3 , and Z3. Other (identical) 
solutions are possible, using the matrix elements yi, zi, and Z 2 - 


4-dimensional case 
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If there is one additional neutrino, decomposition of the neutrino mass diagonalization 
matrix into factors including the PMNS neutrino mixing matrix is more complicated. Lets 
define 2-dimensional rotation matrices in the 4-dimensional complex space, similarly to 
ref. [13], 


n(4) _ 

^12 — 


( Ci2 Si2 0 0\ 
-Si2 Ci2 0 0 
0 0 10 

\ 0 0 0 ly 


p(4) _ 

-^13 ~ 


( Ci3 0 §13 0\ 

0 10 0 

-§13 0 C 13 0 

\ 0 0 0 ly 


n(4) _ 

-^23 ~ 


/l 0 0 o\ 

0 C 23 S 23 0 
0 —S23 C23 0 

\o 0 0 ly 


p(4) _ 

4114 ~ 


/ C 14 0 0 §^4 \ 
0 10 0 

0 0 10 

\ -§14 0 0 C 14 y 


n(4) _ 

2124 ~ 


(l 0 0 0 \ 

0 C 24 0 §24 
0 0 10 
y 0 —§24 0 C 24 y 


p(4) _ 

2134 — 


and phase matrices: = diag fjW _ diag (l, _ 

Note that a shorter notation can be used to define the elements of the rotation matrices: 


/l0 0 0 \ 

0 10 0 

0 0 C34 S^4 

\ 0 0 -§34 C34 y 

r(4) 


(B.ll) 


[^fkVa = + 5 ^ 5 l) + §*, 544 ' - , 


(B.12) 


where 5^ equals 1, when a = 5, or 0, otherwise. This notation is not restricted to the 
4-dimensional case. 

The unitary matrix is parameterized by 


f^(4x4) ^ fy(4) . ^ 




(B.13) 


with the PMNS matrix defined by a product of three rotation matrices: 


L^pmns 0 

0 1 


p(4) p(4) p(4) 
4423 4413 44i2 


(B.14) 


There are 16 parameters in this case, namely: 6 mixing angles (^ 12 , ^ 13 ; 43, ^ 14 ) ^ 4 , 44 ); 
1 Dirac phase hi 3 ; 2 Majorana phases ki and K 2 ] 3 phases Sia, 62 A, 44; and the matrix 
containing 4 phases (j)i = 1) 2,3,4). 

For our model with ur = 1 we calculate the diagonalization matrix Uioop numerically. 
Defining its elements as 


^( 4 x 4 ) 


/ Xi X2 X 3 X4\ 

yi y 2 ys yA 

Zl Z2 Z3 ZA 

\ h t2 t 3 ti y 


(B.15) 
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and comparing to eq. (B.13) we find the relations: 


^*12 

014 


= arcsin 
= arcsin (|x 4 |), 





6*13 = arcsin 



023 = arcsin 



024 = arcsin 
= arg(xi), 



034 = arcsin 
arg(t4), 



(B.16) 

(B.17) 

(B.18) 


(/)2 =arg(xi) - arg(x 2 ) 

aby 2 + X 2 xly4'/b-s/a- |x 3|2 + d|o| X 2 x%^a - |y 4 p/ (a\/c) 


— i In 


^/a -^/hc - \d\‘^ sjh - \x2\‘^ 


(B.19) 


</>3 = </>2 + * In 


d\a\^/ay/bc — \d\^ \/a — IxsPy^c — |2;4p/|(i| 
a?Vhcz^ + aVby/cx 2 , x\ - |y4p + d|a| yl z^^^a - |x 3 P 


(B.20) 


(5i 3 = arg(xi) - arg(x 3 ) - 4>2 - iln 
(5i 4 = </>i - arg(x4), ^24 = </>2 - arg(?/4), S 34 = argiz^) - 4>3, 

where: 

_ 1 I 12 1 . 1 I 12 I 12 

a = I — \X4\ , 0=1 — l^sl — \X4\ , 

c = 1 - |x4p - |?/4p, d = ay 3 + X3X4y4. 



(B.21) 

(B.22) 

(B.23) 


(B.24) 


Because of the discontinuous nature of the square root function in the complex plane, 
yTcy 7 ^ y/x^ in general. Therefore a simplification of the above expressions is limited. 

Due to its origin and the relations between the elements, the expressions do not contain 
all entries of the rotation matrix defined in eq. (B.15). These relations are obtained 

comparing eq. (B.13) with the following matrix elements from eq. (B.15), forming the 
upper-triangular matrix: xi, X2, X3, X4, y2, 2/3, 2 / 4 , ^3, -^ 4 , and 64. Other (identical) 
solutions are possible using the matrix elements zi, Z 2 , D, ^ 2 , and 63 . 


5-dimensional case 

To introduce factorization containing the PMNS neutrino mixing matrix in the 3 + 2 
case, we first define the rotation matrices in the 5-dimensional complex space, similarly to 
ref. [13] 


/ 

C 12 

S 12 

0 

0 

0\ 


/ 

Cl3 

0 

®13 

0 

0\ 


-S 12 

C 12 

0 

0 

0 

II 


0 

1 

0 

0 

0 


0 

0 

1 

0 

0 


-§13 

0 

Cl3 

0 

0 


0 

0 

0 

1 

0 



0 

0 

0 

1 

0 

V 

0 

0 

0 

0 

1/ 


V 

0 

0 

0 

0 

1/ 
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p{5) _ 

^23 ~ 


p{5) _ 

-^24 ~ 


p(5) _ 

^15 ~ 


p{5) _ 

^35 - 


and the phase matrices: 


/I 


0 0 0 o\ 


/ C14 

0 0 

§M 0\ 

0 

C 23 S 23 0 0 

0 ( 5 ) _ 

, -fri4 — 


0 

1 0 

0 0 

0 


S 23 C 23 0 0 


0 

0 1 

0 0 

0 


0 0 10 



§14 

0 0 

C14 0 

VO 


0 0 0 1/ 


V 

0 

0 0 

0 1/ 

/I 


0 0 

0 o\ 


/I 

0 

0 

0 o\ 

0 

C24 0 

§k 0 

II 

0 

1 

0 

0 0 

0 


0 1 

0 0 

0 

0 C 34 

§34 0 

0 


§24 0 

C24 0 


0 

0 - 

§34 

C 34 0 

Vo 


0 0 

0 V 


Vo 

0 

0 

0 1 / 

/ ClE 

0 0 

0§^5V 


/I 

0 

0 

0 0 \ 


0 

1 0 

0 0 


0 

C25 

0 

0 §25 


0 

0 1 

0 0 

II 

0 

0 

1 

0 0 


0 

0 0 

1 0 


0 

0 

0 

1 0 

V- 

§15 0 0 

0 Ci5 j 


Vo 

—§25 0 

0 C25/ 

/I 

0 

0 

0 0 \ 


/I 

0 0 

0 

0 \ 

0 

1 

0 

0 0 

II 

0 

1 0 

0 

0 

0 

0 

C 35 

0 §35 

0 

0 1 

0 

0 

0 

0 

0 

1 0 


0 

0 0 

C 45 §45 

Vo 

0 

—§35 

0 easy 


Vo 

0 0 

—§45 C 45 / 


C/f ^ =diag (^ £*<>3 ^ ^ ^i<P 5 


=diag (l, e*«i/ 2 , 
The unitary matrix is parameterized by 


[,(5x5) ^ [,(5) , ( 45 ) 45 ) 45 )^( 5 ) 


0(5) 0(5) 0(5) 
-^34 -^24 -^14 




with the inclusion of the PMNS matrix 


(B.25) 

(B.26) 


(B.27) 


C^PMNS 0 

0 1 


p(5) p(^) 

-^23 -^13 -^12 


(B.28) 


There are 25 parameters in the 5-dimensional case: 10 mixing angles (012, ^ 13 ; ^ 23 , 
^ 14 ) ^ 24 ) ^ 34 ; ^* 15 , 025)^355^45); 1 Dirac phase Jis; 2 Majorana phases ki and K 2 ; 7 phases 
(^ 14 , ( 524 , ( 534 , (5i5, S 25 , < 535 , < 545 ; and the matrix ' containing 5 phases 4 )1, i = 1,..., 5. Only 
numerical solutions for the parameters are possible. 

A simplification is possible in our analysis. According to the structure of the diagonal- 
isation matrix C/ioop in the 4- or 5-dimensional cases, the 3x3 sub-matrix in the top-left 
corner is dominant. This sub-matrix corresponds to the matrix t/pMNS- In the numerical 
calculation, it suffices to use the expressions B.6 in order to estimate the oscillation angles 
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^ 12 ) ^13 and 023; in the cases of ur = 1 and ur = 2. There is a numerical precision differ¬ 
ence between the approximated angle values and the values obtained using the expressions 

(B.16) in the 4-dimensional case or the numerical solutions in the 5-dimensional case. The 

approximation speeds up the calculations significantly. 
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